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Abstract —This Work presents the application of control with H m problem, together with pole location 
constraint by LMI formulation in a Buck converter. The uncertainties modeling of the system variables uses the 
constraint by a convex polytope. Mathematical software with the implementation of algorithms for LMIs 
solution, as well as software for circuit simulation, were used as an aid for the formulation and validation of the 
control. The result of the robust controller in a transient period is compared to a classic controller, giving 
justification for the quality of the results. 
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I. INTRODUCTION 

Direct current to direct current (DC-DC) converters are 
generally used in many applications due to its ability to 
increase or decrease voltage with high efficiency. There 
are several topologies of these converters, one of them is 
the synchronous buck converter, shown in Fig. 1, which 
only reduces its input voltage. This when compared to 
non-synchronous buck has higher efficiency at full load[l]. 
In order to deliver a near-desired voltage at the output, 
i.e.a voltage that follows a variable reference, it is essential 


to use a controller that keeps the system stable to all 
desired specifications. 

The use of linear matrix inequalities (LMIs) has 
recently been presenting efficient solutions to control 
problems. This is due to the characteristic of 
simultaneously considering various requirements, 
performance constraints and robustness through the 
formulation of LMIs. In this context, research is growing 
in several areas, focusing on control through these 
inequalities [2]. In this work, for example, a robust 
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//oo control is presented, adding constraints on the convex 
region, in which the closed-loop system poles can be 
placed. The main reference is the study presented in [3]. 

To validate this work, two controllers were applied to a 
DC-DC Buck converter (Fig. 1). The/Co control, containing 
the constraints represented by LMIs, is performed by state 
feedback. In this way, two sensors are required, one for 
reading the converter output voltage and one for reading 
the current over the inductor. The second controller uses 
the classic approach of Proportional Integral (PI) control. 
MATLAB/Simulink software was used for system 
simulation. The SeDuMi, YALMIP and ROLMIP 
packages were used in the solution and optimization of the 
LMIs presented in this work. 

The comparison of these two controllers allows us to 
evaluate the advantage that the control with pole region 
constraint considering its characteristic uncertainties, has 
over the classic PI control. Such a comparison is 
interesting for three reasons, the first is that the load 
connected to the converter output may vary within a 
known range in which the converter is designed to work. 
The second is that the capacitance may decrease, e.g. an 
electrolytic capacitor has a liquid dielectric, which may 
degrade over time. The third is that the real components 
have tolerances,e.g.±20%, and could vary according to 
temperature, such as the capacitor and inductor present in 
the converter. 

Finally, the comparison of the results of this work aims 
to present the differences between a controller with 
higher implementation cost, due to the development time 
and hardware cost, and a PI controller that has a less 
expressive cost.This paper is divided in the following 
order. Section II presents the modeling of the Buck 
converter (Fig. 1). Section III introduces the characteristics 
of the //oo controller. Section IV presents the 
implementation of control techniques and the comparison 
of results. Section V presents the conclusions. 

II. SYSTEM MODELING 

Two steps make up the synchronous Buck converter, 
shown in Fig. 1. On the first step of the switching SW1 
conducts and SW 2 shall be open, on the second one SW 2 
conducts and SW 1 stays open. The duty cycle d represents 
the portion of the switching period T s , in which SW1 is 
conducting. The addition of the second load in parallel to 
the resistance R , represented by the source / 0 , represents a 
disturbance in the output current. 

Equating of the inductor voltage to each switching step 
and its weighting of the values by switching parcel results 
in (1) and (2). The equating of the capacitor current to each 


switching step, following the same weighting, results in (3) 
and (4): 


<V L > = ^-.({V g -V 0 -R L .I L ).d.T s 


+ (-V 0 -R L .I L Ul~d).T s ) 


< V, > = L. 


dh 

dt 


V g .d-V 0 -R L .I L . 

+ ( 4 - Y~ Io )- (1 ~ dyTs 


< Ir >= C. 


dVp 

dt 


= h 


R 


■■■■( 2 ) 

....(3) 

....( 4 ) 


From (2) and (4), it is possible to perform alternating 
current(AC) modeling for small signals. In this, the 
average values summed to an AC portion replace each 
variable. The DC terms that represent the permanent 
response cancel each other. Moreover, it is assumed that 
the second order AC terms are small compared to the first 
order ones, to the extent that you can neglectthem [4]. 


Such process is presented in (5) to (13). 

Vg = Vg+V^ .(5) 

Vo = Vo + Vo .(6) 

4 = 4 + 4.(7) 

4 = 4+4 .( 8 ) 

d = d + d . (9) 

d(4+4) _ _ _ „ ^ ^ _ 

L. v L J = V g .d + Vg.d + Vg.d + Vg.d-V 0 -V 0 

-R l T l -R l Sl .( 10 ) 


d(f L ) ___ 

L dt = V 9 d + V 9 d ~ V 0-RL-h .( 11 ) 

dQb + Vo) _ ^ Vo % - ^ 

c.-^-—— = / L + / L - TT - .(12) 

dt _ 1 R R 0 0 v J 

d(Vo) ^ Vo _ 

C.—= I L - -£-/ 0 .( 13 ) 

dt L R u y J 


The model given by (11) and (13) is a type 0, i.e. the 
plant has no integrator. For the steady-state system error to 
be null, for a ramp type input, it is necessary to insert an 
integrator augmented state, shown in (14) and (15). This 


transforms the system into type 1 [5]. 


X 3= f ( V ref ~Vo)-dt . 

. (14) 

d(x 3 ) 


dt ref 0 . 

.( 15 ) 


Considering (11), (13) and (15), it is possible to 
represent them in state space notation, as presented in (16) 
and (17): 
dx(t ) 

—— = A. x(t ) + B u . u(t) + B w . w(t) + B ref . V ref . (16) 


z(t ) = C z .x(t ) + 

Dzu- u(t) + Dzw w(t) .(17) 

, wherein: 
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X(t)= V 0 .(18) 


h 

Vo 

d\ , 


u(t ) = [d] . (19) 

w(t) = [Q] .(20) 

z(t) = [%] . (21) 

Rr 


--F- —r 0 




C 

1 9 . 

\ 

L 

0 
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~rTc 

-l 


o 

oj 


.( 22 ) 


.(23) 
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B re f 


0 

1 

~c 
0 
0 
0 

.11 
= [0 1 
.= [ 0 ].. 
= [ 0 ]... 


.(24) 


0 ]. 


(25) 

.(26) 

.(27) 

.(28) 


D z 
D z 

The disturbance represented by the w vector, from the 
current source 7 0 , represents the behavior of the converter 
output voltage V 0 , for sudden changes in the output current 
of the converter[3]. The parameters V g , L, C e R are 
uncertain and limited to a minimum and maximum value, 
as represented by (29) to (32): 


V g . ,Vg 

^mim Umax] 


P = V g e 

s _ 1 r 1 1 i 

L \-Lmnr I'm in -I 


* = c e 


J max 

1 


-‘min- 

1 


r ’ r . 

max u min 


B 


R.C 


r n 1 r n 
^max- lx max ^mirr ^min 


.(29) 

.(30) 

(31) 

..(32) 


Therefore, matrices A and B u depend on these N = 4 
uncertain parameters, generating a polytope with 2 N = 16 
vertices. Table 1 shows the polytopic representation of 
matrices A and B ni . 


III. ROBUST CONTROL CHARACTERISTICS 

This Section presents the characteristics of the robust 
control developed in this paper. 

3.1 LMIs 

Considering x E as the variable and F t = F/ as 
non-variant matrices. An LMI has the form given in (33): 


m 

Fix) ±F 0 + I F t Xi >0 .(33) 

i=1 

According to [6], there is an LMI feasibility problem 
when a solution x has to be found, such thatF(x) > 0. 

The analysis of LMIs in dynamic systems has over a 
century of study, beginning in 1890, when Lyapunov’s 
theory was published. This theory shows that a differential 
equation, given in (34), is stable, if and only if, there is a 
positive definite matrixP, which satisfies (35): 


^x(t) = Ax{t) .(34) 

A'P + PA < 0.(35) 


Lyapunov proved that for the LMI in (35) to have a 
solution/ 5 > 0, one must choose any matrix Q > Oand 
solve the equation/!'/ 5 + PA + Q = 0. If P > 0 is found, 
there is a solution [2]. 

The sequence of studies on LMIs showed that the 
classic LQR (Linear Quadratic Regulator) problem, 
associated with the Riccati equation, could be represented 
as an LMI problem, through the Schur complement [2]. 
This verification is important because, if possible, several 
problems can be solved together, if represented in LMI 
form, as shown in (36) with k constraints. 

F(x) = diag{F 0 (x),F 1 (x),...,F k (x )} > 0.(36) 

Over time, several algorithms have been developed for 
the solution and optimization of LMIs, such as 
Nemirovskii’s algorithm [2]. Today, with computational 
packages, LMIs have become an efficient tool in many 
engineering areas, including modern control [7]. 

3.2 Quadratic stability 

The existence of a quadratic function, given in (37), 
that satisfies L(x) <0, is a necessary and sufficient 
condition to assume that the linear system, presented in 
(34), is stable. Thus resulting in (38): 

V(x) = x'Px >0, Vx ± 0.(37) 

V(x) = x'(A'P + PA)x <0, Vx ^ 0.(38) 

Then, the system is stable, if and only if, there is a 
positive definite symmetric matrix/ 5 , for whichL(x) < 0. 
To satisfy (38), the conditions in (39) must be met: 

BP > 0; A'P + PA < 0.(39) 

If the uncertainties in A are polytopic, the feasibility of 
the LMI problem must be verified, i.e. the polytope must 
have j vertices. Then, a P > 0solution must be found, such 
that (40) is met: 

AiP + PAi< 0, i = 1,2, ...,j .(40) 

The feasibility of this problem implies that the system 
is stable for every matrix belonging to the polytope. 
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Table 1: Poly topic representation of matrices A and B u 


A, = 


Rl 

Lmin 

1 


u min 
1 




B Ul = 


V a ■ 
Umim 
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Rl 

Lmax 
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u max 
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Rmin-C 
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Rl 
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Rl 
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1 
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0 
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Rl 
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Rl 

Lmax 
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^min-^max 


-1 


0 

0 
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' B u 14 
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Rl 

Lmin 
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^min 
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K max-^min 


B U7 = 


0 

0 
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Rl 

Lmax 
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u max 
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Rmax-^min 


-1 


/ B u 
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Rl 

Lmin 
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Rmax-Cmax 
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u 8 


0 
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1 16 


Rl 

Lmax 
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L max 
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n max-^max 


-1 


, ft 


u 16 


3.3 Quadratic stability for a closed-loop system with 
state feedback 

The LMIs in (39) are adapted for a closed-loop system 
with state feedback u = Kx. Works [3] and [8] simplify 
these LMIs using the following theorem: 

Theorem 3.1 : The closed-loop system with state 
feedback u = Kx is stable, if and only if, there is a 
symmetric matrixl/k E M nxn and a matrixT E M mxn , such 


that: 

W > 0.(41) 

AW 4- WA' 4- B U Y 4- Y'B U ' <0.(42) 

The state feedback gain is obtained by (43): 

K = YW- 1 .(43) 

3.4 Pole placement 


The main motivation for pole clustering, in a specific 
region on the left side of the complex plane, is the transient 
response characteristics of a linear system. A second order 
system, e.g. with poles — (w n ± jw d , is fully characterized 


in terms of the undamped natural frequency w n = |A|, the 
damping rate<( and the damped natural frequencyw d = 
w n s] 1 — ( 2 .By constraining and fixing A in a prescribed 
region, bounds can be put on the characteristic values, thus 
ensuring a satisfactory transient response. The regions of 
interest, that can be defined, include the a-stability regions 
vertical strips, disks, conic sectors, ( Re(s ) < —a), and 
others [9]. 

The prominent region of this project, which has control 
purposes, is characterized by the set S(a,p,6 ), composed 
of complex numbers a+jw d , that satisfy (44) [9]: 
a < —a <0; |a 4- jw d \ < p; tanOo < —\w d \ (44) 
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Fig. 2: Pole region constraint geometry [9]. 


Fig. 2 (a) demonstrates the set S and its variables.The 
comparison of Fig. 2 (a) and Fig. 2 (b) relates certain 
parameters of a control system, such as the minimum 
decay rate a = <w n , the minimum damping ratio ^ = 
cosO , and the maximum undamped natural frequency w d = 
p. sen.Setting the minimum or maximum characteristics of 
each parameterbounds the maximum overshoot, delay 
time, rise time and settling time [9]. 

From [9] a study is made for the representation of the 
constraint, presented in Fig. 2 (a), by means of LMIs. 
Works [3] and [8] simplify these LMIs by the following 
theorem: 

Theorem 3.2 : The closed-loop poles of a system, 
considering the state feedback control law u = Kx, are 
within the region defined by the setS(a,r, 0 ), if there is a 
positive symmetric matrixl/Fand a matrixT, such that: 

AW + WA' + B U Y + Y'B' U + 2 aW <0.(45) 

-pW WA' + Y'B' U 


AW + B U Y 


-pW 


< 0 . 


.(46) 


' sen9(AW + WA' + B U Y + Y’B ’ U ) 
cosd (-AW + WA' — B U Y + Y’B ’ U ) 

cos9(AW-WA' + B U Y-Y’B' U ) 1 
sen9(AW + WA' + B U Y + Y'B' U )\ . 1 ’ 

O ganho de realimenta§ao de estados e dado por (48): 

K = YW~ X .(48) 

3.5 TheWoo norm 

The/Zooiiorm can be represented by the H 2 norm, being 
equivalent to the largest value obtained by this 
one.Therefore, considering lhe(J operatorsbelonging to 
theffoo space, (49) is defined: 

||CJ 2 

IlSiloo = sup xel2 = sup M2=1 ||5J 2 .(49) 

If the operator is a transfer functionG(s), the definition 
in (50) can be used: 

II c II co = sup w o{G(jw )}.(50) 

Considering that<j{G(/w)}is the maximum singular 
value of G(j(o). If G(s) isa SISO (Single Input and Single 
Output) transfer function, the/Co normcorresponds to the 
largest value of the Bode diagram, for example. 

The/Zoo problemconsists of finding a controller, such 
that the/Zoonorm of the closed-loop transfer function G wz 
(wherewis the disturbance andzis the output) is minimized, 
i.e. less than a specified valuey.This value being 
determined using a particular case of the small gain 
theorem. 

The work [2] explains how to calculate the ZZoonorm for 
a closed-loop system with state feedback, determining the 
dual system LMIs, presented by [10] using the following 
theorem: 

Theorem 3.3: A closed-loop system with state 

feedbacks = Kx and-^%- < y, is stable,if and only if, 
ML 

there is a positive defined symmetric matri xW E 
R nxn and a matrixT E R mxn , such that: 

AW + WA' + B U Y + Y'B' U B w WC' Z + Y'D ZU 

B'w -Y 1 D 'zw 

C Z W + D ZU Y D zw —yl 

< 0.. (51) 

The constraint proposed in (41) and (42) is included in 
(51), so Theorem 3.3 ensures quadratic stability. 

IV. SIMULATION RESULTS 

The model obtained in Section II for the Buck 
converter does not take into account its switching 
frequency f s = 100 kHz. In order to prevent the 

control from destabilizing, it is necessary to limit the 
maximum bandwidth of the controller to a frequency 10 
times lower than the switching frequency [3]. 

The purpose of the ZZooControl is to have a maximum 
overshoot of 20% ( %Over ) and a minimum settling time 
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of 10ms ( T e ). The parameters a, 0 e pcan be related to the 
desired system performance through (52) to (55) [5]. 
-log(%0ver / 100) 


? = 


:= 0,45595 .(52) 


^/tt 2 + io^(%0ver/100) 2 

9 = arcos(0 = 62,8739°.(53) 

4 

877.28.(54) 


a 


f-T e 


f<? 


p = 2. 7i. —“ = 62831.85. 
p 10 


(55) 


The synchronous Buck converter always operates 
similar to the continuous conduction mode, since its 
switches do not depend on the current flowing in the 
inductor to polarize, as it does in the non-synchronous 
converter's diode [l].For any R load, there is a peak-to- 
peak output voltage ripple of AV 0 = 12,5mV, and in the 
inductor there is a peak-to-peak current ripple equivalent 
tod/ L = 500mA. 

Such ripples are the result of the switching process, and 
the maximum values are defined in the converter 
design.The controller itself cannot mitigate these ripples. 
For the converter operating at steady state, the minimum 
values of the inductorLand capacitorCare defined 
according to (56) and (57). In these, the duty cycle/) = 
0.5is where the largest ripples of current and voltage occur 
[11]. The inductor Lin practice is not ideal as it has a fixed 
resistance/^. 


V g .D.(l-D) 

j _ y max 
L'min ~ 


Ah-fs 


(56) 


Vg ■ Z). (1 — D) 

u max 


(57) 


mm 8 .L min .AV 0 .f s 2 . 

The values range of the DC-DC Buck converter 
uncertain parameters is shown in Table 2. 

Table 2: DC-DC Buck converter parameters. 


Parameters 

Values Range 

Rl 

o.i n 

L 

[500,800]/!// 

C 

[50,200 ]fiF 

R 

[2,100]/cfi 

v g 

[80,100] V - 


Using MATLAB, the converter system, described in 
(16) and (17), with its polytopic representation of 16 
vertices, shown in Table 1, was implemented for the 
parameter ranges presented in Table 2.With the aid of the 
SeDuMi, YALMIP and ROLMIP packages, the LMIs, 
represented by (41), (45), (46), (47) and (51), were 
optimized.As a result, the state feedback gain, in (58), was 
obtained. The gamma minimization resulted in y = 
4.0809. The block diagram representation, built with 
MATLAB/Simulink, of the Buck converter with state 
feedback control, is shown in Fig. 3. 

K = [0.3389 -0.4435 603.6809].(58) 

The PI controller was designed, for comparison, with 
the parameters presented in Table 3. As a result, the 
proportional and integral gains are shown, respectively, in 
(59) and (60). 


Table 3: PI control design parameters. 


Parameters 

Values 

Margem de Fase 

45.5950° 

fc 

558 Hz 

L 

650/!// 

C 

125 /if 

R 

5 ka 

Vg 

90 V 


K P = 3.5951 x 10" 4 .(59) 

K\ = 1.1827.(60) 


To evaluate the response of the implemented 
controllers, five test scenarios were defined. In the first 
one, the performance of the controllers was compared, 
without considering the uncertainties in the Buck converter 
parameters.In this scenario, the voltage reference was a 
step signal equivalent ioV re f = 60 Fat the initial time, and 
changed ioV re f = 30 V after0.05 5. The output voltages 
obtained for the robust and PI controllers are shown in Fig. 
4. 
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Fig. 3: Synchronous Buck converter with state feedback control. 



0 0.01 0.02 0.03 0 04 0.05 0.06 0.07 0.08 0 09 0.1 

Time (s) 

Fig. 4: Converter output voltage scenario 1. 

The next four test scenarios evaluated the system 
response considering the uncertainties in the system 
parameters.The voltage reference was a step signal 
equivalent to V re f = 60 V at the initial time, and changed 
to V re f = 30 V after 0.08 5. The voltage input was a step 
signal equivalent to V g = 100 V at the initial time, and 
changed to V g = 80 V after 0.12 s. The load R also varied, 
initially it wasR = 100 /cH, changing to R = 
2 HattimeO.04 s , and to R = 10 HattimeO.06 s. These 
settings were used for the test scenarios 2, 3, 4, and 5. 

Considering the values range of the DC-DC Buck 
converter uncertain parameters, shown in Table 2, the 
scenario 2 converter was configured with inductance 
L = 500 uH and capacitance C = 50 uF. Fig. 5 compares 
the scenario 2 control responses. 



Fig. 5: Converter output voltage scenario 2. 


The scenario 3 converter was configured with 
inductance L = 500 uH and capacitance C = 200 uF. Fig. 
6 compares the scenario 3 control responses. 



Fig. 6: Converter output voltage scenario 3. 

The scenario 4 converter was configured with 
inductance L = 800 uH and capacitance C = 50 uF. Fig. 
7 compares the scenario 4 control responses. 
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Fig. 7: Converter output voltage scenario 4. 


The scenario 5 converter was configured with 
inductance L = 800 uH and capacitance C = 200 uF. Fig. 
8 compares the scenario 5 control responses. 



Fig. 8: Converter output voltage scenario 5. 

It can be observed in Figs. 5 to 8 that for all the 
uncertainties considered in the design, the proposed robust 
controller was stable and within the desired performance 
conditions. The robust controller presented satisfactory 
transient responses to the R step loads and the applied 
input voltage V g , in comparison to the classic PI controller 
operating outside its nominal design range. 

V. CONCLUSION 

The set of linear matrix inequalities (LMIs) presented 
in this work allows us to ensure robustness and closed- 
loop poles placement within the desired region, for any 
state space modeled system. In this work, a Synchronous 
Buck converter was modeled and a state feedback gain, 
obtained from the LMIs optimization, was applied. The 
system’s robustness and the desired performance were 
observed through simulation. 

The proposed robust controller achieved a more 
consistent performance over the adopted parameters range, 
in comparison to the classic PI control. As can be observed 
when the output loads are changed, and a step signal is 
applied to the converter’s input voltage, for the minimum 
and maximum values of the adopted inductance and 
capacitance ranges. 
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